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The hghtest Higgs scalar boson mass in supersymmetry can be raised significantly 
by extra vector-like quark and lepton supermultiplets with large Yukawa couplings 
but dominantly electroweak-singlet masses. I consider models of this type that main- 
tain perturbative gauge coupling unification. The impact of the new particles on 
precision electroweak observables is found to be moderate, with the fit to Z-pole 
data as good or better than that of the Standard Model even if the new Yukawa cou- 
plings are as large as their fixed-point values and the extra vector-like quark masses 
are as light as 400 GeV. I study the size of corrections to the lightest Higgs boson 
mass, taking into account the fixed-point behavior of the scalar trilinear couplings. 
I also discuss the decay branchings ratios of the lightest new quarks and leptons and 
general features of the resulting collider signatures. 
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I. 



INTRODUCTION 



The Minimal Supersymmetric Standard Model [1] (MSSM) predicts that the lightest neutral 
Higgs boson, has a mass that can only exceed that of the boson by virtue of radiative correc- 
tions. If the superpartners are not too heavy, then it becomes a challenge to evade the constraints 
on set by CERN LEPII e+e^ collider searches. On the other hand, larger superpartner masses 
tend to require some tuning in order to accommodate the electroweak symmetry breaking scale. 
In recent years this has motivated an exploration of models that extend the MSSM and can raise 
the prediction for m^jO. 

In the MSSM, the largest radiative corrections to m^o come from loop diagrams involving top 
quarks and squarks, and are proportional to the fourth power of the top Yukawa coupling. This 
suggests that one can further raise the Higgs mass by introducing new heavy supermultiplets 
with associated large Yukawa couplings. In recent years there has been renewed interest [2-20] 
in the possibility of a fourth family of quarks and leptons, which can be reconciled with precision 
electroweak constraints with or without supersymmetry. However, within the context of super- 
symmetry, if the new heavy supermultiplets are chiral (e.g. a sequential fourth family), then in 
order to evade discovery at the Fermilab Tevatron pp collider the Yukawa couplings would have to 
be so large that perturbation theory would break down not far above the electroweak scale. This 
would negate the success of apparent gauge coupling unification in the MSSM. Furthermore, the 
corrections to precision electroweak physics would rule out such models without some fine tuning. 

These problems can be avoided if the extra supermultiplets are instead vector-like, as proposed 
in [21-24]. If the scalar members of the new supermultiplets are heavier than the fermions, then 
there is a positive correction to m/jO. As I will show below, the corrections to precision electroweak 
parameters decouple fast enough to render them benign. 

To illustrate the general structure of such models, suppose that the new left-handed chiral 
supermultiplets include an SU{2)l doublet $ with weak hypercharge Y and an SU{2)l singlet (p 
with weak hypercharge —Y — 1/2, and $ and (j) with the opposite gauge quantum numbers. The 
fields ^ and (j) transform as the same representation of SU (3)c (either a singlet, a fundamental, or 
an anti- fundamental) , and $ and 4) transform appropriately as the opposite. The superpotential 
allows the terms: 



where M$ and are vector-like (gauge-singlet) masses, and k and h are Yukawa couplings to the 
weak hypercharge -1-1/2 and —1/2 MSSM Higgs fields if„ and H^, respectively. In the following, 
I will consistently use the letter k for Yukawa couplings of new fields to and h for couplings 
to H^. Products of weak isospin doublet fields implicitly have their SU{2)l indices contracted 
with an antisymmetric tensor e^^ = — e^^ = 1, with the first component of every doublet having 
weak isospin T3 = 1/2 and the second T3 = —1/2. So, for example, = <I>i$2 — *I*2*I'i) with the 
components $1, $2, ^1, and $2 having electric charges Y + l/2, Y — 1/2, —Y + 1/2, and —Y — 1/2 
respectively. 

The scalar members of the new chiral supermultiplets participate in soft supersymmetry break- 
ing Lagrangian terms: 



(1.1) 



Aoft = + b^cfxt) + akHu^4> - ahHa^cpj + c.c. + m||$p + m||^|2. 



(1.2) 
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where I use the same name for each chiral superfield and its scalar component. 

The fermion content of this model consists of two Dirac fermion-anti-fermion pairs with electric 
charges ib(y + l/2) and one Dirac fermion-anti-fermion pair with electric charges ±{Y — 1/2). The 
doubly degenerate squared-mass eigenvalues of the fermions with charge ±(y -|- 1/2) are obtained 
at tree-level by diagonalizing the matrix 



nip 



MfM^ 







(1.3) 



with 



Mf = 



hVd Ma 



(1.4) 



which is assumed to be dominated by the M$ and M<^ entries on the diagonal. Here Vu = v sin (] 
and Vd = V cos P are the vacuum expectation values (VEVs) of the MSSM Higgs fields Hu and H^, 
in a normalization where v ^ 175 GeV. The scalar partners of these have a squared-mass matrix 
given by, in the basis ($, 0, $ ,0 ): 



2 2 , 

777-5 = "^F + 





a-hVd - 

\ QkVu - kfj,*Vd b. 



bl 

hll*Vu ?T7|^ + A_l 

<S 2 ' 2 

i 



a%Vu - kjivd ^ 


777i+ A, 



(1.5) 



where the A^g^g = — qsiv? 9w\cos{2f3)rn\ are electroweak D-terms, with and q the weak 
isospin and electric charge. The scalar particle squared-mass eigenvalues of cq. (1.5) are presumably 
larger than those of their fermionic partners because of the effects of 777|,, 777^, 777^ and 777^, inducing 
a significant positive one-loop correction to 777^0 • If tan /? is not too small, the corrections to 777^0 
are largest if the fe-type Yukawa coupling is as large as possible, i.e. near its infrared quasi-fixed 
point. 

The fermions of charge ±(3^ — 1/2) have squared mass M|, and their scalar partners have a 
squared-mass matrix 



^|M$|2+ 777| + A_l 



Y—- 



—6$ 



— h* 



(1.6) 



These particles do not contribute to 777^^0 except through the small electroweak D-terms, since they 
do not have Yukawa couplings to the neutral Higgs boson. Since that contribution is therefore 
parametrically suppressed, it will be neglected in the following. 

With the phases of and Hd chosen so that their vacuum expectation values (VEVs) are 
real, then in complete generality only three of the new parameters M$, M^, k and h can be 
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simultaneously chosen real and positive by convention. Nevertheless, I will take all four to be real 
and positive below. (I will usually be assuming that the magnitude of at least one of the new 
Yukawa couplings is small, so that the potential CP-violating effects are negligible anyway.) 

In the MSSM, the running gauge couplings extrapolated to very high mass scales appear to 
approximately unify near Q = Munif = 2.4 x 10^^ GeV. In order to maintain this success, it is 
necessary to include additional chiral supermultiplets, besides the ones just mentioned. These 
other fields again do not have Yukawa couplings to the Higgs boson, so their contribution to Am^o 
will be neglected below. 

I will be assuming that the superpotential vector-like mass terms are not much larger than 
the TeV scale. This can be accomplished by whatever mechanism also generates the ji term in 
the MSSM. For example, it may be that the terms M$ and are forbidden at tree-level in the 
renormalizable Lagrangian, and arise from non-renormalizable terms in the superpotential of the 
form: 

after the scalar components of singlet supermultiplets S and S obtain vacuum expectation values of 
order the geometric mean of the Planck and soft supersymmetry-breaking scales. Then M$, < 1 
TeV can be natural, just as for ix in the MSSM. 

In the remainder of this paper, I will discuss aspects of the phenomenology of models of this type, 
concentrating on the particle content and renormalization group running (section 2), corrections to 
m/jO (section 3), precision electroweak corrections (section 4), and branching ratios and signatures 
for the lightest of the new fermions in each model (section 5). 



II. SUPERSYMMETRIC MODELS WITH NEW VECTOR-LIKE FIELDS 

A. Field and particle content 

To construct and describe models, consider the following possible fields defined by their trans- 
formation properties under SU{d)c x SU{2)l x U{1)y'- 

Q = (3,2,l/6), Q = (3,2,-l/6), = (3, 1, 2/3), C7 = (3, 1, -2/3), 

= (3, 1,-1/3), :D= (3,1,1/3), L = (1,2, -1/2), L = (1,2,1/2), 
S= (1,1,-1), :E= (1,1,1), A^= (1,1,0), iV= (1,1,0). (2.1) 

Restricting the new supermultiplets to this list assures that small mixings with the MSSM fields 
can eliminate stable exotic particles which could be disastrous relics from the early universe. In 
this paper, I will reserve the above capital letters for new extra chiral supermultiplets, and use 
lowercase letters for the MSSM quark and lepton supermultiplets: 

= (3,2,1/6), iZi = (3, 1,-2/3), = (3, 1, 1/3), 

A = (1,2,-1/2), e, = (1,1,1), 

iJ„ = (l,2,l/2), i7rf=(l,2,-l/2). (2.2) 
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with i = 1,2,3 denoting the three famihes. So the MSSM superpotential, in the approximation 
that only third-family Yukawa couplings are included, is: 

W = fiHuHd + ytHuq3U3 - ybHdqsds - VrHdhea- (2.3) 

It is well-known that gauge coupling unification is maintained if the new fields taken together 
transform as complete SU{5) multiplets. However, this is not a necessary condition. There are 
three types of models that can successfully maintain perturbative gauge coupling unification with 
the masses of new extra chiral supermultiplets at the TeV scale. 

First, there is a model to be called the "LND model" in this paper, consisting of chiral super- 
multiplets L,L,N,N,D,D, with a superpotential 

W = MlLL + MnNN + MdDD + kNHuLN - hNHoLN. (2.4) 

Here L, L play the role of <1>, $ and N the role of 0, (j) in eqs. (1.1)-(1.6). In most of the following, 
I will consider only the case that the multiplicity of each of these fields is 1, although 1, 2, or 3 copies 
of each would be consistent with perturbative gauge coupling unification. These fields consist of a 
5 -I- 5 of SU{5), plus a pair''^ of singlet fields. The non-MSSM mass eigenstate fermions consist of a 
charged lepton r', a pair of neutral fermions iy[ 2; and a charge —1/3 quark b'. Their superpartners 
are complex scalars f(2, ^^12341 ^'12- The primes are used to distinguish these states from 
those of the usual MSSM that have the same charges. 

Second, one has a model consisting of a 10 -|- 10 of SU{5), to be called the "QUE model" below, 
consisting of fields Q, Q, U, U, E, E with a superpotential 

W = MqQQ + MuUU + MeEE + kuHuQU - huH^QU. (2.5) 

The non-MSSM particles in this case consist of charge +2/3 quarks t'-^ 21 ^ charge —1/3 quark b', 
and a charged lepton r', and their scalar partners i'l 2 3 4; ^'1 2 f( 2- 

Third, one has a "QDEE model" consisting of fields Q,Q, D, D, Ei, Ei {i = 1,2) with a super- 
potential 

W = MqQQ + MuDD + Ms^EiEi + koH^QD - hoHaQD. (2.6) 

Although this particle content does not happen to contain complete multiplets of SU{5), it still 
gives perturbative gauge coupling unification. The non-MSSM particles in this model consist of 
charge —1/3 quarks 6'^ 25 ^ charge +2/3 quark t', and two charged leptons t( 2, and their scalar 
partners 2 3 41 ^12 and t( 2 3 4- 

The field and particle content of these three models is summarized in Table I. 

In reference [24], it is suggested that a model with extra chiral supermultiplets in 5 + 5 + 10 + 10 
of SU{5), or equivalently (if a pair of singlets is added) 16 + 16 of 50(10), will also result in 



Here I choose the minimal model of this type that includes Yukawa couplings of the kind mentioned in the 
Introduction while not violating lepton number. It is also possible to identify the fields N and N, since they are 
gauge singlets, or to eliminate them (and their Yukawa couplings) entirely. 
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Model 


New supcrmultiplets 


New particles 








others 


spin 1/2 


spin 


LND 


L.L 




D.D 




-/ ^/ Jj 

'^l,2,;-i,i '1,2 "l,2 


QUE 


Q,Q 


U,U 


E,E 


A,2 y r' 


'-1,2,3,4 "1,2 "l,2 


QDEE 


Q,Q 


D,D 


-E'l,2,-E'l,2 


h' f' t' 

"1,2 ^1,2 


"1,2,3,4 ''1,2 ^1,2,3,4 



TABLE I: The new chiral supermultiplets and the new particle content of the models discussed in this paper. 
The notation for 0, follows that of the Introduction. 

gauge coupling unification. However, the multi-loop running of gauge couplings actually renders 
them non-perturbativc below the putative unification scale, unless the new particles have masses 
well above the 1 TeV scale. For example, working to three-loop order, if one requires that the 
unified coupling (defined to be the common value of ai and a2 at their meeting point) satisfies 
the perturbativity condition aunif < 0.35, then the average threshold of the new particles must 
exceed 5 TeV if the MSSM particles are treated as having a common threshold at or below 1 TeV 
as suggested by naturalness and the little hierarchy problem. In that case, the new particles will 
certainly decouple from LHC phenomenology. Even if one allows the MSSM soft mass scale to 
be as heavy as the new particles, treating all non-Standard Model particles as having a common 
threshold, I find that this threshold must be at least 2.8 TeV if the new Yukawa couplings vanish 
and at least 2.1 TeV if the new Yukawa couplings are as large as their fixed-point values. While 
such heavy mass spectra are possible, they go directly against the motivation provided by the little 
hierarchy problem. Furthermore, at the scale of apparent unification of ai and a2 in such models, 
the value of 03 is considerably smaller, rendering the apparent unification of gauge couplings at best 
completely accidental, dependent on the whim of out-of-control high-scale threshold corrections. I 
will therefore not consider that model further here, although it could be viable if one accepts the 
loss of perturbative unification and control at high scales. The collider phenomenology should be 
qualitatively similar to that of the LND and QUE models, since the particle content is just the 
union of them. 

B. Renormalization group running 

The unification of running gauge couplings in the MSSM, LND, and QUE models is shown 
in Figure 1. In this graph, 3-loop beta functions are used for the MSSM gauge couplings, and 
nit = 173.1 GcV and tan/3 = 10, and all non-Standard-Model particles are taken to decouple at 
Q = 600 GeV. (The Yukawa couplings kN and Hn in the LND model and ku and hu in the QUE 
model are set to here for simplicity; they do not have a dramatic effect on the results as long as 
they are at or below their fixed-point trajectories.) The running for the QDEE model is not shown, 
because it is very similar to that for the QUE model. Indeed, it will turn out that many features 
of the QUE and QDEE models are similar, insofar as the U + U fields can be interchanged with 
the D + D + E + E fields. This similarity does not extend, however, to the collider phenomenology 
as discussed in section 5. Note that the unification scale, defined as the renormalization scale Q at 
which ai = 0:2, is somewhat higher with the extra chiral supermultiplets in place; in the MSSM, 
Afunif ~ 2.4 X lO^^ GcV, but M^nif 2.65 x 10^^ GeV in the LND model, and M^nif « 8.3 x lO^*' 
GeV in the QUE and QDEE models. The strong coupling misses the unified ai and a2, but by 
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FIG. 1: Gauge coupling unification in the MSSM, 
LND and QUE models. The running is performed 
with 3-loop beta functions, with all particles be- 
yond the Standard Model taken to decouple at Q = 
600 GeV, and mt = 173.1 GeV with tan/? = 10. 



a 




3: 



1 — ' — I — ' — I — ' — I — ' — I — ' — I — ' — r 




6 8 10 12 14 16 
Log^o(Q/GeV) 



6 8 10 12 14 16 
Log^JQ/GeV) 



FIG. 2: Renormalization group trajectories near the fixed point for in the LND model (left panel) and ku 
in the QUE model (right panel), showing the infrared-stable quasi-fixed point behaviors. Here nit — 173.1 
GeV and tan/3 = 10 are assumed. 



a small amount that can be reasonably ascribed to threshold corrections of whatever new physics 
occurs at Munif- 

The largest corrections to 171^0 are obtained when the new Yukawa couplings of the type kN, 
kjj, or ku are as large as possible in the LND, QUE, and QDEE models respectively. These new 
Yukawa couplings have infrared quasi-fixed point behavior, which limits how large they can be at 
the TeV scale while staying consistent with perturbative unification. This is illustrated in Figure 
2, which shows the renormalization group running^ of the ki\f coupling in the LND model and ku 
in the QUE model. The running ol ko in the QDEE model is very similar to the latter (and so 



In this paper, I use 3-loop beta functions for the gauge couplings and gaugino masses, and 2-loop beta functions 
for the Yukawa couplings, soft scalar trilinear couplings, and soft scalar squared masses. These can be obtained 
quite straightforwardly from the general results listed in [25-27], and so are not given explicitly here. 
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FIG. 3: The contours represent the infrared-stable 
quasi-fixed points of the 2-loop renormalization 
group equations in the plane of Yukawa couplings 
{ki,hi) evaluated dX Q — 500 GeV. The allowed 
perturbative regions (defined hy ki, hi < 3 at 
Q = Munif) are to the left and below the contours. 
The long dashed (blue) line corresponds to fc^r, /ijv 
in the LND model. The solid (black) line corre- 
sponds to kij, hjj in the QUE model, and the nearly 
overlapping short dashed (red) line corresponds to 
ko, ho in the QDEE model. Here mt = 173.1 GeV 
and tan /3 = 10 are assumed. The very nearly rect- 
angular shape of these contours reflects the absence 
of direct coupling between the Yukawa couplings in 
the one-loop (3 functions. 



is not shown). In this paper, I will somewhat arbitrarily define the fixed-point trajectories to be 
those for which the extreme Yukawa couplings are equal to^ 3 at the scale Munif where ai and a2 
unify. Then, assuming that only one of the new Yukawa couplings is turned on at a time, and that 
tan /? = 10 with nit = 173.1 GeV and with all new particle thresholds taken to be at Q = 600 GeV, 
the fixed point values also evaluated at Q = 600 GeV are 

LND model: /ctv = 0.765 or Kn = 0.905, (2.7) 
QUE model: % = 1.050 or hu = 1.203, (2.8) 
QDEE model: ko = 1.043 or ho = 1.196. (2.9) 

Turning on both Yukawa couplings at the same time in each model hardly affects the results at all, 
because ki, hi decouple from each other's beta functions at one loop order for each oi i = N, U, D. 
This is illustrated by the very nearly rectangular shape of the fixed-line contours in Figure 3. 

The phenomenology of supersymmetric models is crucially dependent on the ratios of gaugino 
masses. In the MSSM, it is will known that if the gaugino masses unify at Munif, then working 
to one- loop order they obey Mi/ai = M2/a2 = Ms/as = mi/2/aunif, and this relation has only 
moderate corrections from higher-loop contributions to the beta functions. The presence of extra 
matter particles strongly affects this prediction, however. In Table II, the predictions for Mi, 
M2, and M3 at Q = 1 TeV are given for the MSSM, the LND model, the QUE model, and the 
QDEE model. In the latter two cases, I distinguish between the cases of vanishing extra Yukawa 
couplings and the fixed-point trajectories with {ku, hu) = (3, 0) and (ko, ho) = (3, 0), respectively, 
aX Q = Afunif. As before, I have used tan/3 = 10 and mt = 173.1 GeV, and taken all new particle 
thresholds to be at Q = 600 GeV, and assumed for simplicity that the new scalar trilinear couplings 
vanish at Q = Munif. The results will change slightly if these assumptions are modified, but there 



Formally, it turns out that the 2-loop and 3-loop beta functions for these Yukawa couplings have ultraviolet-stable 
fixed points, although these occur at such large values (> 5) that they cannot be trusted to reflect the true 
behavior. Simply requiring the high-scale value of the Yukawa couplings to be somewhat smaller avoids this issue. 




0.2 0.4 0.6 0.8 1 1.2 
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Mi/mi/2 


M2/mi/2 


M3/mi/2 


M2/M1 


M3/M1 


MSSM 


0.41 


0.77 


2.28 


1.88 


5.53 


LND 


0.32 


0.59 


1.75 


1.86 


5.52 


QUE {ku = 0) 


0.097 


0.147 


0.571 


1.52 


5.90 


(%(M,,if) = 3) 


0.109 


0.176 


0.617 


1.61 


5.66 


QDEE {ko = 0) 


0.094 


0.153 


0.572 


1.62 


6.08 


{kDiAUnii) = 3) 


0.107 


0.178 


0.615 


1.66 


5.72 



TABLE II: Gaugino masses Ma/mi/2 for (a = 1,2,3) and ratios of gaugino masses M2/M1 and M3/M1, 
evaluated at Q = 1 TeV in the models described in the text, assuming unified gaugino masses toi/2 at Munif. 



FIG. 4: Running of gaugino masses in the QUE 
model, assuming a unified value mi/2 at M^nii- 
The gluino mass parameter M3 is evolved accord- 
ing to the 1,2, and 3 loop beta functions in the top 
three lines. The 1 and 3 loop beta functions are 
shown for the parameters Mi (bottom two lines) 
and M2, for which the 2-loop and 3- loop results 
are not visually distinguishable. Note the signifi- 
cant running of M3 due to multi-loop efl:ects. 




^2 4 6 8 10 12 14 16 
Log^o(Q/GeV) 

are a couple of striking and robust features to be pointed out about the gaugino masses in these 
models. First, because the unified gauge couplings are so much larger in the models with extra 
matter than in the MSSM, the gaugino masses at the TeV scale are suppressed relative to mi/2 
by a significant amount compared to the MSSM. Secondly, the one-loop prediction for the ratios 
of gaugino masses is very strongly violated by two-loop effects^ in the extended models, which 
were evidently neglected in [24]. For example, in both the QUE and QDEE models, the one-loop 
prediction is that M3 = mi/2, independent of Q, since the one- loop beta function for M3 happens 
to vanish. However, the correct result is that M3 does run significantly, with M^/mii2 reduced by 
some 40% from unity, depending on the Yukawa coupling value. This reflects, in part, the accidental 
vanishing of the one- loop beta function; in contrast, the three- loop contribution to the running is 
quite small compared to the two-loop one. This is illustrated for the QUE model in Figure 4, which 
shows the renormalization-scale dependence of the running gaugino mass parameters Mi, M2, and 
M3 in the QUE model [with {ki/,hi/) = (3,0) at Q = Munif], evolved according to the 1, 2, and 
3-loop beta functions. 

Another notable feature of the extended models is that they permit gaugino mass domination 



Similar effects have been noted long ago in the context of "semi-perturbative unification" [28]. 
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'"u 








MSSM 


2.08 


2.01 


2.00 


0.67 


0.37 


LXD 


1.89 


1.82 


1.81 


0.6.3 


0.35 


QUE {ku = 0) 


1.24 


1.20 


1.19 


0.45 


0.28 


{ku{M^ni{) = 3) 


1.29 


1.24 


1.24 


0.47 


0.30 


QDEE {ko = 0) 


1.24 


1.20 


1.20 


0.45 


0.28 




1.30 


1.25 


1.24 


0.47 


0.30 



TABLE III: Ratios of first- and second- family MSSM squark and slepton mass parameters to mi/2, evaluated 
at Q = 1 TeV, assuming unified gaugino mass dominance at Q = M^nif (wq = Aq =0). 

for the soft supcrsymmctry breaking terms at the unification scale, according to which all soft scalar 
masses and scalar trilinear couplings are assumed negligible compared to the gaugino masses, or 
Aq = 0, rriQ = in the usual mSUGRA language. In the MSSM, this "no-scale" boundary condition 
is problematic if applied strictly, because it predicts that the lightest supersymmetric particle (LSP) 
is not a neutralino. However, in the QUE and QDEE models, the increased size of the gauge 
couplings at high scales gives extra gaugino-mediated renormalization group contributions to the 
scalar squared masses, so that they are safely heavier than the bino-like LSP. For the squarks and 
sleptons of the first two families, this is illustrated in Table III (for the same models as in Table 
II), by giving the ratios of the running masses to the unified gaugino mass parameter mi/2- The 
contributions of the gaugino masses to the new extra squarks and sleptons in the LND, QUE, and 
QDEE models are listed below: 

LND: (mD,m;p,mL,m;^,miv,rn^) = (1.80, 1.80,0.63,0.63,0,0), (2.10) 
QUE: {mQ,mQ,mu,mjj,TnE,m^) = (1.17,1.29,1.25,0.94,0.267,0.299), (2.11) 
QDEE: {mq, m^, mD,m-^, niE, m-^) = (1.30, 1.18, 0.94, 1.24, 0.266, 0.304). (2.12) 

Here I have chosen to display the results for boundary conditions at Munif of kN = hN = and 
for kij = S,hu = and for ko = 3, /iu = 0, respectively. It should be noted that these are all 
running mass parameters, and the physical mass parameters will be different. Also, if there are 
non-zero contributions to the running scalar squared masses and scalar trilinear couplings at Munif , 
the results will of course change. For example, including a non-zero common ttiq, as in mSUGRA, 
will raise all of the scalar squared masses, yielding a more degenerate scalar mass spectrum. 

For the QUE and QDEE models, we see from Tables II and III and eqs. (2.11) and (2.12) 
that the bino mass parameter is well over a factor of 2 smaller than the lightest slepton mass, for 
unified, dominant gaugino masses. Since neutralino mixing only decreases the LSP mass compared 
to the bino mass parameter, the LSP will be a neutralino. In contrast, for the LND models, the 
gaugino mass dominance boundary condition would predict that the scalar component oi N oi N 
(a non-MSSM sneutrino) should be the LSP, and should be nearly massless. In fact, including a 
non-zero Yukawa coupling /i^r or fc^r would give the corresponding scalar a negative squared mass. 
If there is an additional positive contribution to that sneutrino mass, then it can be the LSP, and 
it might be interesting to consider it as a possible dark matter candidate. 

The corrections to the lightest Higgs squared mass considered in the next section depend on 
the scalar trilinear coupling a^j^, a^^, or a^j^ of the type appearing in eq. (1.2). It is therefore 
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Log^o(Q/GeV) Log^o(Q/GeV) 

FIG. 5: Renormalization group running of scalar trilinear couplings Akj^ in the LND model (left panel) 
and Ak^ in the QUE model (right panel), normalized to mi/2, the common gaugino mass parameter at 
the unification scale Munif- The different lines correspond to different boundary conditions at Munif- The 
corresponding Yukawa couplings kpf and kjj are taken to be near their fixed-point trajectories, with kN = 3 
and = 3 at Mu„if. The running of yl^^ in the QDEE model is very similar to that shown here for Ak^. 

useful to note that these couplings have a strongly attractive fixed-point behavior in the infrared 
when the corresponding superpotential couplings k]^, kjj and ko are near their fixed points. To 
illustrate this, consider the quantities 

^fcjv = o-kj^/kN, Ak^j = akjjjkij, Ak^ = ak^/kD, (2-13) 

for the LND, QUE, and QDEE models, respectively. The renormalization group runnings of Akj^ 
and Akjj (each normalized to ?ti-i/2) are shown in Figure 5, for various input values at the unification 
scale. The running of Akj^ in the LND model is seen to have a mild focusing behavior, leading 
to values at the weak scale of —0.1 < Akj^/mii2 < 0.6 for input values at Munif in the range 
~3 < Akf^/mij2 < 3. In the case of Ak^j in the QUE model, one finds an even stronger focusing 
behavior leading to —0.5 < >lfey/mi/2 ^ —0.3 at the weak scale. The running of Ak^^ in the QDEE 
model is very similar (and so is not shown) . It is useful to note that in the cases of Akjj in the QUE 
model and Akj^ in the QDEE model, most of the contribution to the running comes from the gluino 
mass parameter. This will still be true if one does not assume gaugino mass unification, provided 
only that the gluino mass parameter M3 is not very small compared to the bino and wino mass 
parameters Mi and M2. Therefore, the previous results concerning the fixed-point behavior of Akj-, 
and Akij remain approximately valid if mi/2 is replaced by the value of M3 at the unification scale. 

C. Fine-tuning considerations 

One of the primary model-building motivations in recent years is the super symmetric little 
hierarchy problem, which concerns the tuning required to obtain the electroweak scale, given the 
large supersymmetry breaking effects needed to avoid a light Higgs boson that should have been 
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seen at LEP and to evade direct searches for superpartners at LEP and the Tevatron. One way to 

express this problem is to note that the Z boson mass is related to the parameters and vn^n^ 
near the weak scale by: 

1 1 ^ 

- 2"^l = ImI' + + 2^^^^ + 0(1/ tan2 /?), (2.14) 

where AF is the radiative part of the effective potential. Although there can be no such thing 
as an objective measure of fine tuning in parameter space, the cancellation needed between 
and m\r can be taken as an indication of how "difficult" it is to achieve the observed weak scale. 
Large values of —mj^^ require more tuning in this sense. 

In the MSSM, with the gauge and Yukawa couplings taken to be the values of the infamous 
benchmark point SPSla' [29] for a concrete example, one finds 

-mjj^ = 1.82M| - 0.212M| + O.I56M3M2 + O.O23M1M3 - 0.32itM3 - O.OTitMa 

+0.1lA'^ - QMm%^ + 0.36?fi|3 + 0.28m|3 + . . . . (2.15) 

Here 'rnjj^ on the left side is evaluated at the scale Q = 600 GeV, where corrections to AV are 
presumably not too large. The non-MSSM particle thresholds are also taken to be at the same 
scale. The hats on the parameters on the right side denote that they are inputs at the apparent 
unification scale Munif = 2.4 x 10^^ GeV. They consist of gaugino masses Mi,2,3, scalar squared 
masses ifijj^, and Th\^, and At = at/yt- I have neglected to write other contributions with 
small coefficients. Note that the gaugino masses and scalar squared masses are not assumed to be 
unified here. The essence of the super symmetric little hierarchy problem is that after constraints 
from non-observation of the lightest Higgs boson, the charged supersymmetric particles, and from 
the relic abundance of dark matter are taken into account, the remaining parameter space tends 
to yield —rnjj^ 2> rn'^z/'^^ ^^^^ some fine adjustment is needed between —rnjj^ and It was 
noted long ago in ref. [30] that the gluino mass parameter M3 is actually mostly responsible for 
the tuning needed in w-^^, because of its large coefficient as seen in eq. (2.15), and this problem 
can be ameliorated significantly by taking IM3/M2I smaller than unity at Munif. This can easily be 
achieved in non-mSUGRA models. For example, taking IM3/M2I ~ 1/3 produces near cancellation 
between the M| and M| terms with opposite signs in eq. (2.15), yielding a smaller value for rnjj^. 

Now let us compare to the corresponding formulas in the LND, QUE and QDEE models under 
study here. For the QUE model, I find near the fixed point ku = 1.05 with hu = that the most 
significant contributions are approximately: 

-mj^^ = 2.10M| + 0.035M| + 0.019Mf - O.OI4M3M2 - 0.075itM3 - 0.016itM2 
+0.022ife^M3 + 0.014ife^M2 + 0.057 - 0.015itife^ + 0.25i^^ 
-O.Umjj^ + 0.34rn| + 0.27ml^ + 0.47m| + 0.40m|. + . . . . (2.16) 

Again the hats on parameters on the right side denote their status as inpT.it values at Munif, and 
on the left side is evaluated at Q = 600 GeV, which is also where the new particle thresholds 
are placed, and tan/? = 10. The result of eq. (2.16) seems to reflect a worsening of the little 
hierarchy problem, since the contribution to —rnjj proportional to M| is even larger than in the 
MSSM case, while the physical MSSM superpartner masses are actually lower for fixed values of 



13 



the input soft parameters, as can be seen from Tables II and III. This imphes that for a given scale 
of physical superpartner masses, including notably the top squarks that contribute strongly to m^o , 
one will need larger —mij^^, and thus larger and so a more delicate cancellation between the 
two. Note also that since the contribution proportional to M| is positive (and quite small), there 
cannot be a cancellation as in the MSSM for large \M2/M3\. Counteracting these considerations, 
there is the fact that there are large positive corrections to m^o from the new particles, as discussed 
in the following section, so that the top squark masses need not be so large. 

It is interesting to compare with the corresponding result when the new Yukawa coupling kjj is 
instead taken to vanish: 

-mjj^ = 1.14M| - 0.107M| + O.I53M3M2 + O.O22M1M3 - OASGAtMs - 0.090itM2 

+0.125i? - O.TOmjj^ + OMml^ + 0.21ml^ + ... (2.17) 

for ku = 0. Here the impact on fine-tuning is less because the coefficient of M| is reduced, there 
is no large positive contribution from the new scalar soft masses, and the possibility of significant 
cancellation between the gluino and wino mass contributions (if \Ad2/M^\ > 1) is restored. But, 
counteracting this, there is no large positive contribution to m^o from the extra vector-like sector 
when ku = 0. 

Results for the QDEE model are quite similar. At the fixed point with ko = 1.043, I find 

-mjj^ = 2.12M| + 0.034M| + 0.006Mf - O.OI3M3M2 - O.OSSitMa - 0.017itM2 
+0.029ifc^M3 + 0.014ife^M2 + 0.054i^ - 0.027 At Ak^, + 0.12i^^ 
-0.22m^^ + 0.33m| + 0.26ml^ + 0.37m| + 0.39m|) + . . . , (2.18) 

and for k^ = 0, 

-mjj^ = 1.15M| - 0.106M| + O.I54M3M2 + O.O24M1M3 - 0.439itM3 - 0.090itM2 

+0.125i? - O.TOmjj^ + 0.30m|3 + 0.21ml^ + .... (2.19) 

The same general comments therefore apply for the QDEE model as for the QUE model. 
Treating the LND model in the same way, I find for k]\f = 0.765: 

-rriH^ = 1.74M|-0.166M|+0.131M3M2 + 0.020M3Mi-0.33itM3-0.06i^ 
+0.07ifc^M3 + 0.1 lif - OMAtAk^ + 0.05i|^ 

-0.62m|^^ + 0.38m|3 + 0.30m|3 + . . . . (2.20) 

The dependence on the soft parameters in the sector of new extra particles is very slight, due to the 
fact that the fixed-point Yukawa coupling is not too large. Here we see that even at its fixed point 
the LND model is qualitatively quite similar to the MSSM, in that a ratio of IM2/M3I larger than 
1 at the unification scale can reduce —rnjj and therefore mitigate the amount of tuning required 
with \fi\^. For comparison, the result with fcjv = is 



-mjj^ = 1.68M|-0.178M|+0.164M3M^2 + 0.020M^3Mi-0.36itM^3-0.08it^ 
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+0.12if - QMmjj^ + 0.34rfi|3 + 0.26ml^ + ..., (2.21) 

which shows quite similar characteristics. 

Summarizing the preceding discussion, there are two general counteracting effects on the little 
hierarchy problem from introducing vector-like super multiplets with large Yukawa couplings. The 
impact of contributions to —rnjj^ generally tends to worsen the problem, but the additional cor- 
rection to m^o discussed in the next section works to mitigate the problem. (Ref. [24] obtained 
qualitatively similar results, but with quite different numerical details, presumably due to neglect 
of higher-loop contributions to the running of gaugino masses, as noted above.) I will make no 
attempt to further quantify the competition between these two competing and opposite impacts 
on the little hierarchy problem, because there is simply no such thing as an objective measure on 
parameter space, and because there is great latitude in choosing the remaining parameters anyway. 



III. CORRECTIONS TO THE LIGHTEST HIGGS SCALAR BOSON MASS 

The contributions of the new supermultiplets to the lightest Higgs scalar boson mass can be 
computed using the effective potential approximation, which amounts to neglecting non-zero exter- 
nal momentum effects in hP self-energy diagrams. Since m^o is much smaller than any of the new 
particle masses, this approximation is quite good for these contributions. The one-loop contribution 
to the effective potential due to the supermultiplets in eqs. (1.1)-(1.5) is: 

4 

Ay = 27V,^[F(M|) - F(M|^)], (3.1) 
1=1 

where Nc is the number of colors of and and Mp, are the squared-mass eigenvalues of 
eqs. (1.3) and (1.5), and F{x) = x'^[ln{x/Q'^) — 3/2]/647r^. Here Q is the renormalization scale. I 
will assume the decoupling approximation that the neutral Higgs mixing angle is a /? — 7r/2, 
which is valid if m^o ^ 'm\o ■ Then the correction to m^o is 



Amlo = 



1 d 



idvl Vudvu 



cos^ 13 



+ 



dvj Vddvdi ' dvudvd 



+ sin (3 cos /3 - ^ \ AF. (3.2) 



Before presenting some numerical results, it is useful to note a relatively simple analytical 
result that can be obtained if the superpotential vector-like fermion masses are taken to be equal 
(M$ = = Mp) and the soft supersymmetry-breaking non-holomorphic masses are equal {m\ = 
= = = rn?), and the small electroweak D-terms and the holomorphic soft mass terms 
6$ and are neglected. Then, writing 

Mg = Mp + = average scalar mass (3.3) 

X = Ml /Ml (3.4) 

k = ksm.(3, h = hcos(3, (3.5) 

Xk = Ak- n cot /?, Xh = Ah- n tan /3, (3.6) 
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FIG. 6: The functions f{x) and /max(^c) described in the 
text, graphed as a function of ~ the average ratio of 
scalar to fermion masses. 




and expanding to leading order in the normalized Yukawa couplings k and h, one obtains: 



Amlo 



1 



) 



2. Iw 1 
-- 2 - - (1 - - 



xM2^ 3x^ Ux'^M^^ 

Xk{2Xk + Xh)/{^x^M^) 



[I 



1, 



_^(2 _ _ i) _ Xh{2Xh + Xk)/{^xHl^' 

6 X X 



Zx' Ux'^M^ 



(3.7) 



where 



fix) = ln(x) 



1/ Ix/ Ix 



(3. 



It is often a good approximation to keep only the contribution proportional to k'^, corresponding 
to the case where fctan/? ^ h. In that limit, eq. (3.7) agrees with the result given in [24], which 
can be rewritten as simply: 



Am 



47r2 



k^v^ 



sm 



/3 



J_x 

3x' 



12x2M4 



(3.9) 



Note that x is, to first approximation, the ratio of the mean squared masses of the scalars to the 
fermions. A key feature of the result for Am^o is that the contribution of the vector-like particles 
does not decouple with the overall extra particle mass scale, provided that there is a hierarchy x 
maintained between the scalar and fermion squared masses. To get an idea of the impact of this 
hierarchy, the function f{x) is depicted in Figure 6. In the limit of unbroken super symmetry, /(I) = 
0, and f{x) monotonically increases for scalars heavier than fermions [x > 1). The other significant 
feature that could lead to enhanced Am2Q is the mixing parameterized by X^. The maximum 
possible value of the X^ contribution in eq. (3.9) is obtained when X| = 2M2(3x - 1), leading 
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FIG. 7: Estimates for the correctfons to the Higgs mass 
as a function of -y/x, where x = (M^ + m?) /M^ is the ra- 
tio of the mean scalar squared mass to the mean fermion 
squared mass, in the simphfied model framework used in 
eq. (3.9) of the text, using TVc = 3 and k'^v'^ sin" /5 = (190 
GeV)^, corresponding roughly to the QUE or QDEE 
model near the fixed point with reasonably large tan [3. 
The lower line is the no-mixing case = 0, and the up- 
per line is the maximal mixing case Xk = 2M^(?tx — 1). 
The Higgs mass before the correction is taken to be 110 
GeV. 




FIG. 8: Gorrections to m/jo in the QUE model with 
kjj — 1.05, for varying mi^2 with other parameters 
described in the text. Here Mp — 400, 600, and 
800 GeV is the vector-like superpotential fermion 
mass term, and Ms is the geometric mean of the 
new up- type scalar masses. The upper and lower 
lines in each case correspond to Ak = — 0.5mi/2 
and Ak ~ —0.3ini/2, respectively, at the TeV scale. 
The value of m^jo before these corrections is as- 
sumed to be 110 GeV. 
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to a "maximal mixing" result Am^o = ^^^f ^ sin^ /3/max(2;) where /max = f{x) -|- (3 — l/x)^/3. 
This function is also graphed in Figure 6 to show the maximal effects of mixing from the new 
fermion sector. In Figure 7, I show an estimate of the corresponding corrections to Am/jO, taking 
Nc = 3 and /c*^u^sin^/3 = (190 GeV)^ (corresponding roughly to the QUE or QDEE model near 
the fixed point with reasonably large tan /?) and assuming that the predicted Higgs mass before 
the correction is 110 GeV, so that Am/^o = y^(110 GeV)^ + Am^g - 110 GeV. 

The previous depiction may be too simplistic, since the superpotential and soft supersymmetry- 
breaking masses need not have the simple degeneracies that were assumed. Also, as found in 
the previous section, the scalar trilinear coupling has a fixed point behavior that implies that 
the mixing is neither maximal nor zero (but closer to the latter). A more realistic estimate 
is therefore as depicted in Figure 8. Here, I take scalar masses inspired by the renormaliza- 
tion group solutions of the previous section for the QUE model. In particular, I take three 
cases for the vector-like superpotential masses at the TeV scale, Mq = Mjj = Mp = 400, 600, 
and 800 GeV. The Yukawa couplings are taken to be at the fixed-point values ku = 1.05 and 
hu = ^- The soft supersymmetry-breaking terms are parameterized by {mQ,mQ,mu,m-jj) = 
(1.17, 1.29, 1.25, 0.94)mi/2 and (6q,6c/) = -{Mq, Mu)mi/2, and Ak = -0.3mi/2 and -0.5mi/2, 
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all at a rcnormalization scale of 1 TeV. The results turn out to be not very sensitive to bg or 
bjj, or to the MSSM supersymmetric Higgs mass parameter fi (taken to be 800 GeV here), or to 
tan/3 as long as it is not too small (tan/3 = 10 was used here). Figure 8 shows the results for 
Am/jO = ^ (110 GeV)^ + Am^o ~ HO GeV, as a function of Ms, the geometric mean of the scalar 
masses. Quite similar results obtain for the QDEE model at the fixed point with ko = 1.043, 

hD=0. 

Figure 8 illustrates that the contribution of the new extra particles to m^o is probably much less 
than the "maximal mixing" scenario, if one assumes that the TeV-scale parameters (particularly the 
scalar trilinear coupling Uk) can be obtained by rcnormalization group running from the unification 
scale. Note that the models illustrated in Figure 8 represent gaugino-mass dominated examples. 
If one assumes that the soft scalar squared masses at Munif actually have significant positive 
values (from, for example, running between Munif and Mpianck), then the low-scale model will be 
even closer to the no-mixing scenario, since the diagonal entries in the scalar mass matrix will be 
enhanced, while the mixing terms are still subject to the strong focusing behavior seen in Figure 
5. 

In the case of the minimal LND model, one expects the maximum contributions to Am^o 
to be suppressed by a factor of roughly {kN/kuf/N,- ^ (0.765/1.05)^/3 ^ 0.094. This leads 
to corrections that are typically not large compared to the inherent uncertainties in the total 
prediction. This counts against the minimal LND model as a way of significantly increasing the 
Higgs mass. One can also consider n> 1 copies of the LND model, with each k^ Yukawa coupling 
near a common fixed point to maximize Am^o- However, then the common fixed point value is even 
smaller, with fcjv = 0.695 for n = 2, and fcjv = 0.650 for n = 3. (A much more significant correction 
to m^o can occur if one enhances the model with several copies of the extra fields connected by 
the "lateral" gauge group idea of [23].) 



IV. PRECISION ELECTROWEAK EFFECTS 

Because the Yukawa couplings responsible for large effects on m^o break the custodial symme- 
try of the Higgs sector, it is necessary to consider the possibility of constraints due to precision 
electroweak observables arising from virtual corrections to electroweak vector boson self-energies. 
In this section, I will show that these corrections are actually benign (and much smaller than 
previously estimated), at least if one uses only Mj, My^/, and Z-peak observables as in the LEP 
Electroweak Working Group analyses [31, 32] rather than including also low-energy observables as 
in [33]. The essential reason for this is that the corrections decouple with larger vector-like masses, 
even if the Yukawa couplings are large and soft supersymmetry breaking effects produce a large 
scalar-fermion hierarchy. Indeed, they decouple even when the corrections to m^o do not. 

The most important new physics contributions to precision electroweak observables can be 
summarized in terms of the Peskin-Takeuchi S and T parameters [34]. For the measurements of 
Standard Model observables, I use the updated values: 

sle = 0.23153 ± 0.00016 ref. [31] (4.1) 

Mw = 80.399 ± 0.025 GeV ref. [32, 35] (4.2) 

= 83.985 ± 0.086 MeV ref. [31] (4.3) 

a4^^(Mz) = 0.02758 ±0.00035 ref. [31] (4.4) 
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Mt = 173.1 ± 1.3 GeV ref. [36] (4.5) 
as{Mz) = 0.1187 ±0.0020 ref. [33] (4.6) 

with Mz = 91.1875 GeV held fixed. For the Standard Model predictions for s^g, Mw, and in 
terms of the other parameters, I use refs. [37], [38], and [39], respectively. These values are then 
used to determine the best experimental fit values and the 68% and 95% confidence level (CL) 
ellipses for S and T, relative to a Standard Model template with Mt = 173.1 GeV and = 115 
GeV, using 

1 + -^S-^T, (4.7) 



(■5eff)sM 4:S^C2W C2W 

- ^- " + (4.8) 



-— — = 1 - adwS + a(l + sl^diy)^, (4.9) 
(1 ^ jSM 

where sw = sin^^y, cw = cos0w, S2W = sm{26w), C2w = cos{26w), and dw = (1 — 4sp^^)/[(l — 
4s^ + 8s^)c2w]- The best fit turns out to be 5 = 0.057 and T = 0.080. 

The new physics contributions to S and T are given in terms of one-loop corrections to the 
electroweak vector boson self-energies II^vi^, II^^, 11^^ and 11^,^, which are computed for each 
of the LND, QUE and QDEE models in Appendix A. They are dominated by the contributions 
from the fermions when the soft super symmetry-breaking scalar masses are large. It is useful and 
instructive to consider the simplified example that occurs when, in the notation of the Introduction, 
M$ = = Mp with an expansion in small = kvu, rrid = hvd, and Mw- Then one finds for 
the new fermion contributions: 

= [4iml + ml) + m„md(3 + 20Y^)] , (4.11) 

where is the weak hypercharge of the left-handed fermion doublet, denoted $ in the Introduction, 
that has a Yukawa coupling to (so that Y$ = —1/2, 1/6, and —1/6 for the LND, QUE, and 
QDEE models respectively). Equations (4.10), (4.11) agrcc^ with the results found in [40, 41]. 
An important feature of this is that the corrections decouple quadratically with increasing Mp, 
regardless of the soft supersymmetry breaking terms. This is in contrast to the contributions 
to Am^o, which do not decouple as long as there is a hierarchy between the scalar and fermion 
masses within a heavy supermultiplet. It also contrasts with the situation for chiral fermions (as 
in a sequential fourth family), which yields much larger AS, AT. 



^ However, note that the result for AT quoted in ref. [24] actually corresponds to the improbable case hvd = kvu, 
rather than h = 0. So, for small h, the actual correction to AT is almost a factor of 4 smaller than their estimate. 
As a result, much smaller values for Mp are admissible than would be indicated by ref. [24]. 
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FIG. 9: Corrections to electroweak precision observ- 
ables S, T from the LND model at the fixed point 
ikN,hN) = (0.765,0), for varying Ml = Mn = 
rrir' > 100 GeV, in the hmit of heavy scalar super- 
partners. The seven dots on the line segment corre- 
spond to rrir' = 100,120,150,200,250,400 GeV and 
00, from top to bottom. The experimental best fit is 
shown as the x at {AS, AT) = (0.057,0.080). Also 
shown are the 68% and 95% CL elhpses, obtained 
as described in the text. The point AS* = AT = 
is defined to be the Standard Model prediction for 
mt = 173.1 GeV and m?,o = 115 GeV. 
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If h = 0, then the results of eqs. (4.10), (4.11) from the fermions become, numerically: 

AT = 0.54Ar, sm\p) (^^^^|^) ' , (4.12) 
AS = O.lSiV, sm\p) ( ^°°^^"^ ) ' • (4.13) 

These rough formulas show that it is not too hard to obtain agreement with the precision elec- 
troweak data, provided that Mp is not too small, but it should be noted that especially for light 
new fermions with mass of order 100 GeV, the expansion in large Mp is not very accurate, with 
eqs. (4.12) and (4.13) overestimating the actual corrections. 

A more precise evaluation, using the formulas of Appendix A, is shown in figures 9 and 10, 
which compares the experimental best fit and 68% and 95% CL ellipses to the predictions from 
the models. Note that in these figures I do not include the contributions from the ordinary MSSM 
superpartners, which are typically not very large and which become small quadratically with large 
soft supersymmetry breaking masses. Figure 9 shows the corrections for the LND model at the 
Yukawa coupling fixed point {kj\f, h^) = (0.765,0), for varying M^r = Mp = m^' > 100 GeV as a 
line segment with dots at m,-' = 100,120,150,200,250,400 GeV and oo. These contributions are 
due to the fermions u'-^2^'^'^ with their scalar superpartners assumed heavy enough to decouple. 
Note that in the LND model b' and b'^ 2 do not contribute to S, T as defined above, since they 
do not have Yukawa couplings to the Higgs sector. Figure 9 shows that even for ttIt-/ as small as 
100 GeV, the S and T parameters remain within the 68% CL ellipse, and can even give a slightly 
better fit to the experimental results provided that mr' > 120 GeV. If the Yukawa coupling kN is 
less than the fixed point value, or if Ml < Mm, then the corrections to S and T are smaller, for a 
given m^/. 

Figure 10 shows the corrections for the QUE model at the Yukawa coupling fixed point 
{ku,hu) = (1.050,0), for varying M{/ = Mq = niy as a line segment with dots at m^/^ = 
275, 300, 350, 400, 500, 700, 1000 GeV and 00. [For a comparison to the approximate formulas (4.12) 
and (4.13), the appropriate values are Mp = my ~ 355, 381, 432, 483, 584, 786, 1088 GeV and 
00, respectively.] Here, I have included the contributions from the scalar states ^1,2,3,4 and 61^2, 
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FIG. 10: Corrections to electroweak precision observ- 
ables S, T from the QUE model at the fixed point 
iku,hu) = (1.050,0), for varying Mq = Mu ^ mi,,, 
with mi/2 — 600 GeV and Ak^, — — 0.4mi/2 and 
mo = and bg ^ bjj = — rni/2MQ, using eqs. (2.11) 
and (2.13). The eight dots on the line segment corre- 
spond to TOt/ = 275, 300, 350, 400, 500, 700, 1000 GeV 
and oo, from top to bottom. The experimental best 
fit is shown as the x at (AS", AT) = (0.057,0.080). 
Also shown are the 68% and 95% CL ellipses, ob- 
tained as described in the text. The point AS* = 
AT = is defined to be the Standard Model pre- 
diction for nit — 173.1 GeV and rUho = 115 GeV. 
Results for the QDEE model are very similar. 
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obtained for mi/2 = 600 GeV and = — 0.4mi/2 and mo = and bq = bu = —mi/2MQ, using 
eqs. (2.11) and (2.13). Smaller values of would imply a chargino lighter than the LEP2 bound; 
see Table II. From Figure 10 we see that a slightly better fit than the Standard Model can be 
obtained for m^/^ > 400 GeV, but even for m^/^ as light as 275 GeV, the corrections remain within 
the 95% CL ellipse. The corrections to S and T for a given m^/ are even smaller (and so the fit 
is even better) if any of the following conditions apply: the Yukawa coupling kjj is below its fixed 
point value, mi/2 or uiq is larger so that the new squarks are heavier, or Mq ^ Mjj. For t'l masses 
less than about 400 GeV, the fit also improves slightly if m^o is larger than 115 GeV. 

I have also looked at the QDEE model at its fixed point {ko,h£)) = (1.043,0), with scalar 
squared-mass soft terms given by eq. (2.12) with = 600 GeV. The results are nearly identical 
to those found in the QUE model in figure 10 with the values for m^/^ replaced by my^, and so are 
not depicted. 

If one also included lower energy data as used in [33], the fits to S and T would be somewhat 
worse, so it is important to keep in mind that the results above are sensitive to the choice of 
following the LEP Electroweak Working Group [31, 32] in using fits based on the Z-pole data and 
mt and m\Y. With this caveat, one may conclude that the models considered here fit at least as 
well as the Standard Model, provided that the new quarks with large Yukawa couplings are heavier 
than roughly 400 GeV, and can do even better if the new squarks are heavy enough to decouple. 



V. COLLIDER PHENOMENOLOGY OF THE EXTRA FERMIONS 

The extra particles in the models discussed above will add considerable richness to the already 
complicated LHC phenomenology of the MSSM. A full discussion of the different signals, and how 
to disentangle them, is beyond the scope of the present paper, but it is likely that the most impor- 
tant distinguishing collider signals will arise from production of the new fermions, especially the 
new quarks. This is simply because of the relatively large production cross-section compared to 
the scalars, which are presumably much heavier due to the effects of soft supersymmetry break- 
ing masses. One can therefore expect signals from direction pair production of the lightest new 
quark, and possibly also from cascade decays of somewhat heavier fermions down to them. For 
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concreteness, I will concentrate on only the final states from decays of the lightest new quark in 
each model. In general, the lightest new quark and the lightest new lepton would be stable, were 
it not for mixing with the Standard Model fermions. At least some small such mixing is necessary 
to avoid a cosmological disaster from unwanted heavy relics. If the mixing is very small, then the 
new fermions could be quasi-stable, with decay lengths on the scale of collider detectors. Then the 
collider signatures will involve particles that leave highly-ionizing and slow tracks in the detectors, 
or feature macroscopic decay kinks or charge-changing tracks. These can be either the new charged 
leptons or hadronic bound states of the new quarks. Such signals have been discussed before in a 
variety of different model-building contexts; for some reviews, see refs. [49-52]. 

In the following, I will assume that the mixing of the new fermions with Standard Model fermions 
is large enough to provide for prompt decays. Mixing of the new fermions with the first and second 
family Standard Model fermions is highly constrained by flavor-changing neutral currents, since the 
vector-like gauge quantum number assignments eliminate the GIM-type suppression. Therefore, I 
will assume that the mixing is with the third Standard Model family, for which the constraints are 
much easier to satisfy. Then the final states of the decays will always involve a single third-family 
quark or lepton, together with a W, Z, or boson. Below, I will discuss the possibilities for the 
branching ratios of the new quarks and leptons, and their dependence on the type of mixing. 

There are existing limits on the extra quarks coming from Tevatron, although these have mostly 
been found with assumed 100% branching ratios for particular decay modes (which as we will see 
below is not necessarily likely). The current limits are, for prompt decays: 

• mt' > 311 GeV for BR{t' Wq) = 1, based on 2.8 fb^^ [42] 

• mb> > 325 GeV for BR(6' ^ Wt) = 1, based on 2.7 fb-^ [43] 

• rub' > 268 GeV for BR(6' ^ Zb) = 1, based on 1.06 fb'^ [44] 

• mb' > 295 GeV for BR(6' ^ Wt,Zb,h^b) = 0.5,0.25,0.25, based on 1.2 fb"^ [45] 
and for quasi-stable quarks: 

• mt' > 220 GeV, based on dE/dx for 90 pb"^ at ^/s = 1.8 TeV [46] 

• rub' > 190 GeV, based on dE/dx for 90 pb"^ at ^/s = 1.8 TeV [46] 

• nib' > 170 GeV for 3mm < cw < 20mm, based on 163 pb~^ [47] 

Also, if the cross-section upper bound found from time-of-flight measurements with 1.0 fb^^ in 
ref. [48] for stable top squarks also applies to stable t' quarks with no change in efficiency, then I 
estimate a bound mt' > 360 GeV should be obtainable, with a somewhat weaker bound for stable 
b' due to a lower detector efficiency. 

At hadron colliders, the production cross-section of the new quarks is due to gg and qq initial 
states and is mediated by the strong interactions, and so is nearly model-independent when ex- 
pressed as a function of the mass. The leading order cross-section is shown in figure 11 for the 
Tevatron pp collider at ^/s = 1.96 TeV and for the LHC pp collider with -y/s = 7, 10, 12, and 14 
TeV. Note the Tevatron will probably be unable to strengthen the existing constraints very signif- 
icantly, at least for promptly decaying new quarks, due to the rather steep fall of the production 
cross-section with mass. At the LHC pair production of mostly vector-like quarks should provide 
a robust signal; see for example studies (in diverse other model contexts) in refs. [6, 54-58]. (Note 
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FIG. 11: Production cross-section for new 
quarks as a function of the mass, for the Tcva- 
tron pp collisions at y/s = 1.96 TeV, and for 
the LHC pp collisions with y/s — 7, 10, 12, 
and 14 TeV. The graph was made at lead- 
ing order using CTEQ5L0 parton distribu- 
tion functions [53] with Q = niqi and apply- 
ing a K factor of 1.5 for LHC and 1.25 for 
Tevatron. 
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that in the models under study here, there is no reason that the flavor-violating charged current 
couplings should be large enough to enable a viable signal from single q' production in association 
with a Standard Model fermion through t-channel W exchange, unlike in other model contexts as 
studied in refs. [59-65].) The branching ratios and possible signals for the LND, QUE, and QDEE 
models are examined below. 



The LND model 



In the LND model, the fermions consist of a b', r', and two neutral fermions u'^ and • The v'^ 
is always lighter than the r'. The fermions b' and v'^ can therefore decay only through their mixing 
with the Standard Model fermions from the superpotential 



(5.1) 



where e^, e^, and are new Yukawa couplings that are assumed here to be small enough to 
provide mass mixings that can be treated as perturbations compared to the other entries in the 
mass matrices. 

First consider the decays of b'. The mass matrix for the down-type quarks resulting from 
eqs. (2.3), (2.4), and (5.1) is: 



^oVd VbVd 



(5.2) 



with eigenstates b and b' . The b' decay can take place only through the coupling, to final states 
Wt, Zb, and hPb. Formulas for these decay widths are given in Appendix B. To leading order, the 
branching ratios only depend on the mass of the 6', and the results are graphed in Figure 12. Note 
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m^, [GeV] m^, [GeV] 

FIG. 12: The branching ratios of the hghtest new quark b' (left panel) and the lightest new lepton ly'i (right 
panel) in the LND model. The !/[ results assume that £„ ^ e^; if instead 3> then BR{v[ — > Wt) — 1 
(not shown). 



that in the limit of large m?,/ , the branching ratios are "democratic" between charged and neutral 
currents, approaching 0.5, 0.25, and 0.25 for Wt, Zh, and h^h respectively, in accord with the 
Goldstone boson equivalence theorem. However, for smaller masses, kinematic suppression reduces 
the Wt branching ratio, so that, for example, the three final states have comparable branching 
ratios for my in the vicinity of 300 to 400 GeV. 

The LHC signals include pp b[b'i W~^W~ti ^ W'^W^W^W^bb. When two same-charge 
Ws decay leptonically and the other two W^s decay hadronically, this leads to a same-charge 
dilepton plus multi-jets (including two b jets) plus missing transverse energy signal, with a total 
branching ratio as high as 25%. This signal is also the basis for the current Tevatron bound 
nib' > 325 GeV, but this assumes BR(6' — > Wt) = 100%; since the actual branching ratio predicted 
by the LND model for that mass range is more than a factor of 3 smaller, the model prediction 
for the signal in the channel that was searched is more than an order of magnitude smaller, and 
decreases sharply for lower rrif,/^. In over half of the other b[b'i production events, there will be four 
or more b jets, coming mostly from events with h^b bbb decays but also from Zb bbb. The 
Tevatron limit [44] of rribi > 268 GeV from assuming BR(6' Zb) = 100% is in a mass range where 
the actual branching ratio is about 0.55, so the actual predicted signal from the LND model is 
more than a factor of 3 smaller. The limit of nib' > 295 GeV from [45], a search which is motivated 
in part by [66, 67], is based on the idealized large mass limit "democratic" branching, but in the 
relevant mass range the model prediction has BR(6' Wt) more than a factor of 2 smaller, and 
decreasing very rapidly for smaller rub', due to the kinematic suppression. The neutral current 
decays, including Z could also play an important role at the LHC, see for example [58] 

for a similar case. 

The decay of z/J in the LND model is dependent on two different mixing Yukawa couplings 
ejv and e^. The mass matrix for the neutral leptons in the {L, N, i^, L, N) basis resulting from 
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eqs. (2.3), (2.4) and (5.1) is 







Ml 



where M^, 



^ Ml hNVd ^ 



kNVu Mm CnVu 



(5.3) 



with the masses of the Standard Model neutrinos neglected. The corresponding mass eigenstates 
are a Standard Model neutrino u and two extra massive neutrino states u[ and fg. The mass 
matrix for the charged leptons is 



Me 



^ -Ml esVd^ 



V 



VrVd 



(5.4) 



Formulas for the resulting decay widths for iy[ Wt and Zu and h^v are given in Appendix B. 
If one assumes that S> ejv, then the decay u[ Wt has a nearly 100% branching ratio. If the 
opposite limit applies, 3> e^, then the branching ratios as a function of mj,/ are as shown in 
the right panel of Figure 12. Note that in the limit of large m^'_^, the branching ratios for Wt, 
Zv, and Ip^v asymptote to 0.5, 0.25, 0.25 respectively when e^v dominates, again in accordance 
with Goldstone boson equivalence with equal charged and neutral currents. So, depending on 
which Yukawa coupling dominates, one could have interesting hadron collider signatures from 
v{v{ production, such as W^W't^t' , and h^h^ + Ef^, and ZZ + Ef^, and Wh^ + Ef^ 
and Z]r^ + E^^^^. So far, there are no published limits specifically on nii^/ based on collider pair 
production with these final states. If < Mjv in this model, then r' will be not much heavier 
than v'l, and so there will be additional contributions to the signal from t'v'i production and r'+r'" 
production, followed by r' —>■ W^*'>i''. It should also be noted that production of u'^ 2 and t[ might 
well be dominated by cascade decays from heavier strongly interacting superpartners. 



B. The QUE model 

In the QUE model, the lightest of the new quarks is always the charge 2/3 quark t'l- After being 
pair-produced at hadron colliders, it can decay due to mixing with the Standard Model fermions 
through the superpotential 

W = euH^qsU + e'^H^Qu^ - eoHaQd^, (5.5) 

where e^, e'^, and Cd are new Yukawa couplings that are assumed here to be small enough to treat 
as perturbations compared to other entries in the mass matrices. The resulting mass matrices for 
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FIG. 13: Branching ratios for the hghtest extra quark, t[, in the QUE model with Mq — Mij, to final states 
Wb, Zt, and h'^t, as a function of m^^ . The left panel shows the "democratic" case that arises when 
dominates (with equal charged and neutral currents), and the right panel shows the 'W-phobic" (mostly 
neutral current) case that arises when e'^ dominates. In the 'W-philic" case that arises when eo dominates, 
then BR{t' Wb) = 1 (not shown). 



the up-type quarks and down-type quarks are 



^ Mq kuVu e'ui^y 

huvd Mu 
y euVu ytVu ) 



^ -Mq eoVd^ 



VhVd 



(5.6) 



with mass eigenstates t,t[,t2 and b,b' respectively. Formulas for the resulting decay widths for t[ 
to Wb, Zt, and h^t are presented in Appendix B. I will concentrate on the three cases where one 
of the mixing Yukawa couplings in eq. (5.5) dominates over the other two. The branching ratios 
depend on the mass of t'^ and on the type of mixing. If e^, provides the dominant effect, then 
the decays are dominantly charged-current, or 'W-philic", with BR(t']^ Wb) = 1. This is the 
scenario for which the Tevatron limit is now > 311 GeV [42]. If instead e'jj dominates, then 
the decays are dominantly neutral-current, or 'W-phobic"; in the limit of large m^/, the branching 
ratios asymptote to BR{t[ Wb) = and BR{t[ Zt) = BR{t[ h^t) = 0.5. Finally, if 
dominates, the the decays are "democratic", with branching ratios for Wb, Zt, and approaching 
0.5, 0.25, and 0.25 respectively in the large m^/ limit. Numerical results are shown in Figure 13 
as a function of m^i^, for the case that ku is at its fixed point value, and hjj = 0, and Mq = Mjj. 
(The results are only mildly sensitive to the last two assumptions.) By taking the different mixing 
Yukawa couplings e^, e'y, and e^, to be comparable, one can get essentially any result one wants for 
the branching ratios, but it seems reasonable to assume that one of the individual mixing Yukawa 
couplings dominates in the absence of some organizing principle. So the possible signatures will 
include W^W~bb, (similar to the Standard Model ti signature, but with larger invariant masses; 
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FIG. 14: Branching ratios for t' decays to Wiy, 
Zt, and h'^r in the QUE and QDEE models, as a 
function of to,-'. 
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see [6, 54, 55, 57, 58] for recent studies of comparable signals), and ZZti and h^h^tt, etc. If 
Mq < Mjj in this model, then the b' will be not much heavier than the t', and one should expect 
an additional component of the signal from b'F and b'b' production, followed by b' W^*H'. 

The r' in the QUE model mixes with the Standard Model r lepton through a superpotential 
term: 



W = -e^HdhE. 



(5.7) 



The mass matrix for the charged leptons resulting from this and eqs. (2.3) and (2.5) is: 



Me 



^ Me ^ 



(5.8) 



with mass eigenstates r and r'. It follows that r' can decay to Wu, Zt, and h^r, with decay 
widths that are computed in Appendix B. Because there is only one relevant Yukawa mixing term, 
the branching ratios depend only on m,-/. They are shown in Figure 14, assuming m/jO = 115 GeV. 
The largest branching ratio for r' is always to Wi', and in the large iBt-/ limit, Goldstone boson 
equivalence provides that the Wu, Zt, and h^r branching ratios approach 0.5, 0.25, and 0.25, 
respectively. The most immediately relevant searches at hadron colliders will be in the mass range 
of rriT-i just above 100 GeV, where the electroweak pair-production cross-section can be sufficiently 
large, and limits do not presently exist. However, note that the appearance of r{ could easily be 
dominated by cascade decays from heavier strongly interacting superpartners. 



C. The QDEE model 

In the QDEE model, the new fermions consist of a b[, b2, t', and r{, Tg. In this model, the 
lighter charge —1/3 quark b'^ is always lighter than the t'. The decays of b'^ in the QDEE model 
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are brought about by superpotential mixing terms with third-family quarks: 



W 



(5.9) 



In the gauge eigenstate basis, the resulting mass matrices for the down-type quarks and up-type 
quarks are 



Md 



^ Mq kDVu ^ 

hovd Md e^Vd 
ye'jjVd ybVdJ 



Mu 



^ -Mq ^ 



(-uVu ytVu 



(5.10) 



with mass eigenstates 6,6^,62 and t,t' respectively. Formulas for the resulting decay widths for b'l 
to Wt, Zb, and h^b are given in Appendix B. As in the case of the QUE model, I will consider 
the three cases where one of the mixing Yukawa couplings in eq. (5.9) dominates over the other 
two. Then the branching ratios depend on the mass of b[ and on the type of mixing. If 
provides the dominant effect, then the decays are dominantly charged-current, or 'W-philic", with 
BR(6'^ Wt) = 1 provided that it is kinematically allowed. The resulting signal at hadron colliders 
will be b'lb'i W~^W~ti W^W^W~W~bb. This is the scenario for which the Tevatron limit 
is presently m^' > 325 GeV [43], based on the same-charge dilepton plus 6-jets signal already 
mentioned above for the LND model. If instead is dominant, then the decays are dominantly 
neutral-current, or 'W-phobic", with BR(6'^ — Wt) = 0; in the limit of large rub', the branching 
ratios slowly approach BR(&'^ Zb) = BR(6'^ — h^b) = 0.5, but with h^b larger for finite masses. 
Finally, if e^, is dominant, the the decays are "democratic", with branching ratios for Wt, Zb, and 
h^b approaching 0.5, 0.25, and 0.25 respectively in the large nib'^ limit. The predicted branching 
ratios are shown in Figure 15 as a function of my^ for the latter two cases, assumings is at its 
fixed point value, and hn = and Mq = Md- (However, it should be noted that, unlike in the 
QUE model case, the results shown are somewhat sensitive to the last of these assumptions.) Note 
that in the "democratic" case, the branching ratios are similar to what one obtains for the b' of 
the LND model. The CDF limit rrib' > 295 GeV was obtained in the idealized case of branching 
ratios obtained in the high mass limit, but for finite m^/, the actual BR(6' — ^ Wt) is much smaller 
and BR(6' — h'^b) is larger. In contrast, the same-charge dilepton signal from b'ib[ — W'^W'~tt is 
turned off in the "Vl^-phobic" case, where the largest overall branching ratio is typically to six b 
quarks, yielding the interesting possible signal b'^b'^ h?h?bb bbbbbb. Decays to Icptons through 
Z bosons are unfortunately suppressed by both small BR(Z i~^i^) and small BR(6']^ —> Zb) 
in this case. If Mq < Md in this model, then the t' will be not much heavier than the 6', and 
one should expect an additional component of the signal from t'b' and t't' production, followed by 
t' ^ vf(*)6'. 

For t[ in the QDEE model, the branching ratio situation is essentially the same as for the QUE 
model as discussed above. 
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FIG. 15: Branching ratios for the hghtest extra quark, h'^^ in the QDEE model with Mq = Md, to final 
states Wt, Zb, and h°b. The left panel shows the "democratic" case that Ed dominates, and the right panel 
shows the 'W-phobic" case that e'^ dominates, leading to mostly neutral-current decays. In the "VF-philic" 
case that eu dominates leading to mostly charged-current decays, then BR{b[ Wt) = 1 (not shown). 



VI. OUTLOOK 

In this paper, I have studied super symmetric models that have vector-like fermions that are 
consistent with perturbative gauge coupling unification and have large Yukawa couplings that can 
significantly raise the Higgs mass in supersymmetry. Some of the more important features found 
for these models are: 

• There are three types of models consistent with perturbative gauge coupling unification and 
all new particles near the TeV scale. The first type (LND) contains up to three copies of 
the 5 + 5 of SU{5). The second type (QUE) contains a 10 + 10 of SU{5). The third type 
(QDEE) is not classifiable in terms of complete representations of SU (5), but consists of the 
fields Q, D, E, E and their conjugates. 

• A complete vector-like family (i.e. a 16 -|- 16 of 50(10)) could also be entertained, but 
was not considered here because a multi-loop renormalization group analysis shows that this 
would forfeit perturbative unification and high-scale control unless (at least some of) the 
new particles are much heavier than 1 TeV. 

• The constraints imposed by oblique corrections to electroweak observables are rather mild, 
especially in comparison to the corresponding constraints on a chiral fourth family, and are 
easily accommodated by present data as long as the new quarks with Yukawa couplings are 
heavier than about 400 GeV, and perhaps considerably lower. 

• The model framework is consistent with the hypothesis that gaugino masses dominate soft 
supersymmetry breaking near the unification scale, without problems from sleptons being 
too light as is the case in so-called mSUGRA models. 
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• The lightest Higgs mass can be substantially raised in the QUE and QDEE models if the 
Yukawa couplings are near their fixed points. However, the extent of this is limited if one 
takes seriously the prediction for the fixed point behavior of the scalar trilinear couplings, 
which limits the mixing in the new squark sector. For example, if the new quarks are 
at Mp = 400 GeV, and their scalar partners have an average mass of Ms = 1000 GeV, 
then one finds an increase in m/jo of up to about 15 GeV (see Figure 8). For larger Ms, 
this contribution increases, but at the expense of apparently more severe fine-tuning of the 
electroweak scale. 

• Despite the sizable positive contribution to the lightest Higgs, the contributions to the jj. 
parameter are also raised, so it is difficult to make any unambiguous claim for an improvement 
in the supersymmetric little hierarchy problem. 

• The new fcrmions can decay through any mixture of neutral and charged currents to third- 
family fermions and W, Z, hP weak bosons, but with different combinations correlated to the 
possible superpotential couplings that mix the new fermions with the Standard Model ones. 

• Existing bounds from direct searches at the Tevatron do not significantly constrain the pa- 
rameter space of these models after precision electroweak constraints are taken into account. 

The collider phenomenology of the MSSM augmented by the new particles in these models should 
be both rich and confusing, leading to a difficult challenge at the LHC and beyond in deciphering 
the new discoveries. 



Appendix A: Contributions to precision electroweak parameters 

This Appendix gives formulas for the contributions of the new chiral supermultiplets to the 
Peskin-Takeuchi precision electroweak parameters [34] . For convenience I will follow the notations 
and conventions of [68] . The oblique parameters S and T are defined in terms of electroweak vector 
boson self-energies by 



aS 



qT 



Hzz(M|)-Hzz(0) 



C2W 

Cy/Sw 
2 



nz^(M|)-n^^(M|) /Ml 



= nw(o)/M^ - nzz(o)/M|. 



(A.l) 
(A.2) 



In the following, the one-loop integral functions G{x), H{x,y), B(x,y), and F(x,y) are as defined 
in ref. [68] , and particle names should be understood to stand for the squared mass when used as 
an argument of one of these functions, which also have an implicit argument s which is identified 
with the invariant mass of the self-energy function in which they appear. 

1. Corrections to electroweak vector boson self-energies in the LND model 

For the LND model, define the gauge eigenstate new neutral lepton mass matrix by 



Ml knVu 



hNVd Mn 



(A.3) 
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and unitary mixing matrices L and R by 



L*MuR^ = diag(mj,/,mj,/), 



(A.4) 



and note m,-/ = Ml- Then the {i^\,i^2,t') fermion contributions to the electroweak vector boson 
self-energies are: 



N, 



167r2 
167r2 



167r2 

2 



(|^,,tP + b5.,tP)G(r') 



(A.5) 
(A.6) 

(A.7) 



An 



WW 



■ ■ -I - J -J " J •' J 

16^2 E [(l<r'tl' + l4.'tl')^(r^^0-4Re(4,,t4,,t)m.'m,,S(r^^.0], (A.8) 

i=l,2 



where Nc = 1 and Ce = — 1 and the massive vector boson couphngs with the new leptons are: 



2cw 



fr/,-/t 



2cw 



RtiRji, 



Z Z 9 f ^ 2 



(A.9) 

(A.IO) 
(A.ll) 



To obtain the (t'i^2,3,4' ^1,2) scalar contribution, consider the new sneutrino squared-mass matrix: 



m| + Aio bl al^Vu-nkNVd^ 



Ml = Ml + 



m 



N 



bL O'hN'^d - IJ'*hNVu rri^ + ^- 1 ,0 

\akj^Vu - l^kNVd bN 



'N 



2 



, (A.12) 



N 



where the supersymmetric part (also equal to the fermion squared-mass matrix) is: 



( t \ 



mIMu 



(A.13) 
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and miy = Mq. Then the {t[,t2,b') fermion contributions to the electroweak vector boson self- 
energies are: 



An 



■77 



AUz^ 
AUzz 



An 



WW 



167r2 
167r2 



gsw 



el E Git'i) + elG{b') 

i=l,2 



167r2 

2 



j=l,2 
,Z i2 I |„Z |2\ 



12 Z Z 



(A.24) 
(A.25) 



167r2 



(A.26) 



E [(l^tT^'tl' + \9^in\')Hib',t'^ - iRe{g^,„g^^„)m,,m,Bib',t-)] , (A.27) 



i=l,2 



where Nc = 3 and = 2/3 and ea = —1/3 and the massive vector boson couphngs with the new 
quarks are: 



4'.t 

^ J 



9b'b't 

ft'6't 



cw \2 

z _ 9 



1 2 

edSw 



9pft — — 



1 



(A.28) 

(A.29) 
(A.30) 



To obtain the {t[ 2 3 4) ^'1 2) scalar contribution, consider the up- type squark squared-mass ma- 
trix: 



M? = + 

u u ' 



m 



Ai 2 

2 ' 3 







ahu^d- IJ'*huVu rn^ + -2 







bu 







, (A.31) 



where the supersymmetric part (also equal to the fermion squared-mass matrix) is: 



^MuMi ^ 



MlMu 



(A.32) 
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Also, the down-type squark mass matrix is 



a 



MA + m^ + A_i_i 



V 



M^ + m^ + Ai 1 

V Q 2 '3 



(A.33) 



Now define unitary scalar mixing matrices U and V by 



UMIu^ = diag(m| , m| , m| , m| ) , VM^V^ = diag(m? , m? ) . (A.34) 



"1 "2 "3 "4 



Then the scalar contributions to the vector boson self-energies are: 

4 2 



An 



■77 



An 



Z7 



Anz2 = 



An 



167r2 
167r2 

167r2 
AT, 



c 2 2 
5 Sw 



gsw 



eiY,ni'iA)+eiY^F{yrh\) 

1=1 1=1 

4 2 



i=l 



i6.2EEk^*i^^(^^>*~^)> 
1=1 j=i 



(A.35) 
(A.36) 

(A.37) 

(A.38) 



where the vector boson couplings with the new squarks are: 



z _ _9_ 

Vj cw 



gZv, = 9{V*iUji-V^U,s)/V2. 



91' b'* = ^d^wj ^ij, (A.39) 

(A.40) 



3. Corrections to electroweak vector boson self-energies in the QDEE model 
For the QDEE model, define the gauge eigenstate new down-type quark mass matrix by: 



and unitary mixing matrices L and R by: 



(A.41) 



R*MdL^ = diag(mb;,m^). 



(A.42) 
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and note m^/ = Mq. Then the (6^,62)*') fermion contributions to the electroweak vector boson 
self-energies are: 



An 



■77 



AUz^ 
AUzz 



167r2 
167r2 



el E + elG{t') 



i=l,2 



gsw 



i=l,2 



(A.43) 
(A.44) 



167r2 

2 



(ki,tr + i5|ptr)G(t') 



An 



WW 



167r2 



(A.45) 



E [d^^Tt'tl' + l55'j;.tP)^(^',fc^) - iRe{g^,,,g^pt)mt,m,,Bit',b'^] , (A.46) 



i=l,2 



where Nc = 3 and = 2/3 and = —1/3 and the massive vector boson couphngs with the new 
quarks are: 



1 3 



z 
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cw 



9m 
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Cw K"^ 



z 
'9t't'f 



cw V2 



-9L*i/V2, 



.w 



9i'.i'\ - 5-^21/ V 2. 



]^RaRji + CdS^Si^ , (A.47) 

(A.48) 
(A.49) 



To obtain the (fe'i 2 3 4) '^'i 2) scalar contribution, start with the down- type squark squared- mass 
matrix: 



M| = Mj + 



Ai 1 
2 ' 3 





bQ 



b*Q al^Vu - fJ-koVd 

aho'^d - fJ'*hDVu rriq + 

^2 



ho 



, (A.50) 



where the supersymmetric part (also equal to the fermion squared-mass matrix) is: 



Mi 
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MdMl 
^ M^Md j 



(A.51) 



35 



Also, the up-type squark squared-mass matrix is given by: 



-bQ 



Now define unitary scalar mixing matrices U and V by: 



UMjU^ = diag(m? , m| , m? , m| ) , VM~V^ = diag(m| , m| ; 

Then the scalar contributions to the vector boson self-energies are: 
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An, 
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where the vector boson couplings with the new squarks are 



z 

9t'*b'. 



9 

cw 

g{V*,Uji - V*,U,s)/V2. 
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cw 



(A.52) 



(A.53) 



(A.54) 
(A.55) 

(A.56) 

(A.57) 



-^ + euSw^ Sij, (A.58) 
(A.59) 



Appendix B: Formulas for decay widths of new quarks and leptons 

This Appendix gives formulas for the decay widths of the lightest of the new quarks and leptons 
to Standard Model states. These decays are assumed to be mediated by Yukawa couplings that 
provide small mass mixings that can be treated as perturbations compared to the other entries in 
the mass matrices. In the following, A(x, y, z) = x'^ + y'^ + z'^ — 2xy — 2xz — 2yz. 

1. Decays ofh' in the LND model 

In the LND model, the lightest quark b' can decay to Standard Model states because of the 
mixing Yukawa parameter in eq. (5.1). In terms of the mass matrix Md in eq. (5.2), define 
unitary mixing matrices L and R by: 



L*MdR^ = diag(mb,mb/). 



(B.l) 
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The relevant couplings of h' to Standard Model particles are 

w T* / /o z 9 



9b'bt 



2cw 



22-^12) 



yb'b = -sm{a){ybRi2 + eDRii)L22/V2, 
- sm{a){ybR22 + eoi?2i)i^i2/v^- 



%6 



(B.2) 

(B.3) 
(B.4) 



It follows that the decay widths of b' are: 

nb'^Wt) = ^\g^^^\^xy\l,rw,rt){l + n-2rw + {l-rtf/rw), 



(B.5) 
(B.6) 
(B.7) 



I z ,2, 

where mb is neglected for kinematic purposes and ri = mf/ m^, for i = Z,W,hP. 

2. Decays of v[ in the LND model 

Consider the decays of the lighter new neutral lepton in the LND model, brought about by 
the superpotential mixing terms e^v and in cq. (5.1). Define unitary mixing matrices L (3 x 3) 
and (2 X 2) in terms of the neutral lepton mass matrix in eq. (5.3) by: 



/ 



R*MlL^ = 



mjy/ 



0m 



(B.^ 



A) 



where we are neglecting the tau neutrino mass. Also define unitary matrices L' and R' in terms of 
the charged lepton mass matrix in eq. (5.4) by: 



L'*MeR'^ = diag(mr, 



Then the relevant couplings of ^ to Standard Model particles are: 

9'^^ri = 9{L*2,L[, + L*2sL[2)/V2 = gRl,R!i,/V2 



z 

9u[u'< 



9 



■(L21L11 + ^23-^13) 



2cw 

IjO cos a , ^ , ^ , „ sin a , ^ „ 
Vu'u = — y^{eNLi3 + kNLu)Ri2 y=-hNLi2Rn. 



V2 

It follows that the decay widths of u[ are: 



V2 



(B.9) 

(B.IO) 
(B.ll) 

(B.12) 



T{u[^Wr) = ■-^{l-rwf{2 + l/rw){\9^^,A' + \9^^,A% 



TTl ' 



(B.13) 
(B.14) 
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m,/ 



rK ^ h'ur) = ^(1 - rnomlj', (B.15) 

where rrir and m^.^ are neglected for kinematic purposes and = mf/m^, for i = Z,W,h^. 

3. Decays of t\ in the QUE model 

Consider the decays of t'^, the hghtest new quark in the QUE model, brought about by the 
superpotential mixing terms in eq. (5.5). Define unitary mixing matrices L, R, L', R' in terms of 
the mass matrices in eq. (5.6) by: 

L*MuR^ = diag(mt,m(/^,m(^), L'*MdR'^ = diag(m6,m6/). (B.16) 

Then the relevant couplings of t[ to Standard Model particles are: 



5ti6t = ff(^2iAi + i23^i2)/V2, gl^, = gR*,,R'n/V2, (B.17) 

^ 2^^^^'^^^^'^^^^^^^^' ^ ~2^"^^^"^^^' (B.18) 

Q COS cx sin ex 

= — (ec/-L23-Rl2 + eyL21-Rl3 + %i^21-Rl2 + 2/t^23-Rl3) -^/lC/-L22-Rll, (B.19) 

= {euLi3R22 + e'^Liii?23 + kuLnR22 + ytLiiR2i) - ^%Li2-R2i. (B.20) 

It follows that the decay widths of are: 

^ = " ^^^'^^ + lA^)(l5jf6tl' + bgtP), (B.21) 

r(tl ^ Zt) = ^aV2(i,^^,^,) [(1 + r,- 2tz + (1 - n)Vrz)(l4,tl' + b|ftl') 

+12VF^Re(5|,t5|tt)]' (B-22) 
= '^>^^'\^^^h^^n)[{l + n-rt,.)Ul^^ + \y'(,\^^ (B.23) 

where the bottom quark is treated as massless for purposes of kinematics and = m|/m^/ for 
i = t,Z,W,hP. 

4- Decays of b'l in the QDEE model 

Consider the decays of 6'^, the lightest new quark in the QDEE model, brought about by the 
superpotential mixing terms in eq. (5.9). Define unitary mixing matrices i?, L, i?', L' in terms of 
the mass matrices in eq. (5.10) by: 

R*MdL^ = diag{mb,mb'^,mi,0, R'*MuL'^ = diag(mt, m*/). (B.24) 

Then the relevant couplings of 6'^ to Standard Model particles are: 

5?;t = 9{L*2iL'n + L*M/V2, g^^p = gRl,R[,/V2, (B.25) 
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9b{b^ = -^^(-^21-^11 + -^23-^13), ^K't'^ ~ 2^^^^^^' (B.26) 

Q sin Q^ cos cx 

Vy^l = — ^ {eoL2iRi2 + e'j,L2iRn + hDL2iRi2 + 2/6^23-^13) + —^kDL22Ru, (B.27) 

Q Sin (y. cos cy. 

y\f, = {eDLizR22 + e'j,LnR23 + hDLnR22 + ybLizR2z) + —^kDLi2R2i- (B.28) 

It follows that the decay widths of b[ are: 

r(6i ^ Wt) = ^aV2(i, rw, rt) [(1 + n- 2rw + (1 - nf /rw){\9^^,, \' + \g^j, f) 

+12V^e(<7feYtt55[ft)], (B.29) 
^ = " ^^^'^^ + V^z)(l46tl' + \9§^lA% (B.30) 

where the bottom quark is treated as massless for purposes of kinematics and rj = m|/m^, for 
i = t,Z,W,h^. 

5. Decays of t' in the Q UE and QDEE models 

Consider the decays of r' in the QUE model, brought about by the superpotential mixing term 
€e in eq. (5.7). In terms of the mass matrix eq. (5.8), define unitary mixing matrices L and R by: 

L*MeR^ = diag(mT-, m-r'). (B.32) 

Then the relevant couplings of r' to Standard Model particles are: 

ffXt = 5^22/^^, af'rl = -^^22^12, (B.33) 

y^"^ = -smia)L22iyrRi2 + eERii)/V2, (B.34) 
y^!^ = -sm{a)Li2iyrR22 + eER2i)/V2. (B.35) 

It follows that the decay widths of r' are: 

T{r'^Wu) = ^[l-rw)\2 + l/rw)\g^j', (B.36) 
r(r' ^ Zr) = 1^(1 - rzn2 + l/rz)\9^,^,\' , (B.37) 
r(r' - h'r) = 1^(1 - r,o f{\y!;:U' + ly^U'), (B.38) 

where = mf/m'^i for i = Z,W,h^, and m,- is neglected for kinematic purposes. In the QDEE 
model, the same calculation holds, provided that Me is replaced by M^-^ corresponding to the 
lighter mass eigenstate nir'- 
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Note added: shortly after the present paper, one with some related subject matter appeared [69]. 
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